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^ Abstract 
*— > 

O | It has recently been suggested that black holes may be described as con- 

densates of weakly interacting gravitons at a critical point, exhibiting strong 
quantum effects. In this paper, we study a model system of attractive bosons 
in one spatial dimension which is known to undergo a quantum phase transi- 
tion. We demonstrate explicitly that indeed quantum effects are important at 
the critical point, even if the number of particles is macroscopic. Most promi- 
nently, we evaluate the entropy of entanglement between different momentum 
modes and observe it to become maximal at the critical point. Furthermore, 
we explicitly see that the leading entanglement is between long wavelength 
modes and is hence a feature independent of ultraviolet physics. If applicable 
to black holes, our findings substantiate the conjectured breakdown of semi- 
classical physics even for large black holes. This can resolve long standing 
^ mysteries, such as the information paradox and the no-hair theorem. 

'X 

1 Introduction 
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In a recent series of papers [H 121 El H] , Dvali and Gomez have proposed a new con- 
ceptual framework for the understanding of black hole physics. The first of two 
key claims is that the black hole is a bound state or condensate of many weakly 
interacting (i.e. long-wavelength) gravitons. Secondly, it was suggested that this 
condensate is at a quantum critical point and therefore exhibits properties that are 
not apparent in the traditional description in terms of (semi-) classical general rela- 
tivity. Most importantly, the underlying quantum physics could be able to resolve 
the mysteries of the information paradox. Hawking evaporation is described as the 
depletion and evaporation of the condensate and its purification is thus a natural 
result. To the same extent, black holes could carry quantum hair [3]. These effects 
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are not visible in the semiclassical approximation, since this limit corresponds to an 
infinite number of black hole constituents. 

In this letter, we will review the Dvali-Gomez proposal and will elaborate in 
more detail why one can expect Bose condensates at a critical point to display 
qualitatively new phenomena. In particular, we will discuss how quantum physics 
can be relevant on macroscopic scales in such systems. To this end, we are going 
to investigate in detail the quantum phase transition of the attractive Bose gas in 
1 + 1 dimensions. 

The idea that black holes may be described by the quantum physics of N weakly 
interacting gravitons was first put forward in [I]. There it was observed that in 
pure Einstein gravity a black hole of mass M and hence size R s = G^M primarily 
consists of gravitons of wavelength A ~ R s . As each of these long wavelength 
particles contributes energy E% = h/X, one obtains 

Ar MM 2 

N ~ — = . (1) 

E l Ml 1 ] 

The idea of a black hole as a Bose condensate of gravitons can also be motivated 
in a bottom-up approach. As gravitons are self coupled, they can potentially form 
a self sustained bound state. The properties of such a bound state can be estimated 
via the virial theorem, 

(£kin> ~ (V) . (2) 

The kinetic energy E^i n of N gravitons of wavelength A is given by 

<£kin> = , (3) 

while a naive estimate of the potential energy of the configuration of size R is 

(4) 

Assuming the size to be of the order of the wavelength, R ~ A, one obtains 

A ~ VNL P . (5) 

It is easily verified that this relation is nothing but Eq. 0. An order of magnitude 
estimate of graviton-emission gives a result consistent with the rate of emission of 
Hawking radiation. Consequently, a self sustained bound state of gravitons, if it 
exists, will likely behave like a black hole. 

One should observe that the typical interaction strength between two gravitons 
is a ~ l/(XMpi) 2 . However, all mutual interactions add up and their total effect 
should be quantified by Not. This implies that the self bound graviton condensate 
is at aN ~ 1, where interactions start to dominate over the kinetic term. This 
condition characterizes the critical point of a zero temperature phase transition 
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or quantum phase transition (QPT) in a simple bosonic model system j3] and is 
considered to ensure that quantum effects are important even for macroscopic black 
holes. 

To substantiate this idea, it is our interest to gain more qualitative insight into 
bosonic systems at a critical point by a detailed study of this 1 + 1-dimensional 
nonrelativistic attractive Bose gas on a ring. The transition in this system was 
discovered and first studied in [5j . We will substantiate the existence of this critical 
point by studying appropriate characteristics. 

We will then focus on the quantum behavior. As a measure of quantumness, 
we calculate the entanglement of different momentum modes applying analytical as 
well as numerical techniques. We observe that it becomes maximal at the critical 
point and for low momentum modes. We interpret this as further evidence that 
the black hole condensate picture can be successful independent of the ultraviolet 
physics that completes Einstein theory. 

The remainder of the paper will be organized as follows. In section 2 we will 
introduce in detail the 1 + 1-dimensional attractive Bose gas, remind the reader of 
mean g and introduce the basis of our numerical studies. Further evidence for the 
existence of a quantum critical point is provided in section 3. We will then introduce 
the fluctuation entanglement as a relevant measure of quantumness and present our 
results in 4. Finally, in the conclusions, we discuss the qualitative consequences of 
our findings with regards to the physics of macroscopic black holes. 



2 The 1 + 1-dimensional Bose Gas 

Throughout this paper, we consider a Bose gas on a ID-circle of radius R with 
attractive interactions. The Hamiltonian is given by 



1 f 2n 



(6) 



where a is a dimensionless, positive coupling constant. This Hamiltonian can be 
cast into a more convenient form by decomposing ip{9) in terms of annihilation 
operators: 

k=— oo 

which leads to 

oo oo 

H = } j k 2 a\a k - — a\a\a m+k ai- m (8) 

k=— oo k,l,m=—oo 

Note that in order to improve readability we have now switched to units R = h = 
2m = 1. The total number operator is 

N= / d9ft(9)4>(6)= 4^- ( 9 ) 
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It was first shown in [3] that an increase of the effective coupling aN on the ring 
leads to a transition from a homogenous ground state to a solitonic phase, where 
the critical point is reached for aN = 1. 



2.1 Mean Field Analysis 

A mean field approach to the hamiltonian ^ leads to the Gross-Pitaevskii energy 
functional 

E[* GP ] = ^ d9 \\d e ^(6)\ 2 - ^ W)| 4 1 (10) 
Jo 1 l \ 

The ground state wavefunction \l/o is obtained through minimization of the energy 
functional subject to the constraint J d6\ty(6)\ 2 = N. This leads to the time inde- 
pendent Gross-Pitaevskii equation 

[d 2 + 7ra|^o(#)| 2 ] Mo) = t*M0) > (ii) 

where /i = dE/dN is the chemical potential. Solutions to this equation are given by 
(see e.g. f6jf] 

M9) ^{S^ d «(^ e -^ m ) ^^^ M ' (12) 

Here, 9q denotes the center of the soliton and m is determined by the equation 

7T\2 



K(m)E(m) = ( - ) aN . (13) 



For small aN < 1, (10) is minimized by the homogenous wavefunction. On the 
other hand, for aN > 1 the solitonic solution has a lower energy. At aN = 1, both 
configurations are degenerate in energy - a clear indication for a quantum phase 
transition. 

On a side note, one may wonder whether the one-soliton solution is stable for 
arbitrary aN > 1 or if multi-soliton solutions may eventually be energetically fa- 
vored. This can be checked in a simple argument. A soliton of size R s has a total 
energy 

N N 2 

£~^-a-. (14) 

Minimization with respect to R yields R s = and E\ = —^a 2 N 3 . A split 
into two stable solitons of boson number rN and (1 — r)N yields a total energy 
Ei = — ja 2 iV 3 [l — 3r(l — r)]. This is bigger than E\ for any r < 1. This can 



^^Here, dn(u|m) is a Jacobi elliptic function and K(m) and E(m) are the complete elliptic 
integrals of the first and second kind, respectively. 
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be straightforwardly generalized two multi-soliton solutions; therefore, the single 
soliton is stable. 

Finally, let us note that the apparent spontaneous breaking of translation sym- 
metry in the solitonic phase is in no contradiction to known theorems about the 
absence of finite volume symmetry breaking. The Gross-Pitaevskii ground state 
only becomes exact in the iV — > oo limit. In this limit, translated Gross-Pitaevskii 
states are orthogonal and do not mix under time evolution. Technically, symmetry 
breaking is made possible because expecation values of composite operators made 
out of the fields diverge in the large N limit. We comment on this in more detail in 
the Appendix. 

This again emphasizes how the classical limit really emerges as a large N limit 
from quantum mechanics. Exactly how this argument breaks down at the critical 
point and what the implications of this breakdown are will be the focus of the 
remainder of this manuscript. 

2.2 Bogoliubov Approximation 

The Gross-Pitaevskii equation is the zeroth-order equation in an expansion of the 
field operator into its mean value and quantum (and, in more general setups, ther- 
mal) fluctuations around it: 

$(9) = $(9)) + 6${6) . (15) 

The spectrum of these small excitations around the mean field can then be found 
in the Bogoliubov approximation. Generally, this corresponds to approximating 
the fluctuation Hamiltonian by its quadratic term and subsequent diagonalization 
through canonical transformations of the field. 

For a>N < 1, i.e. on the homogeneous background, it is convenient to stick to 
the momentum decomposition ^ and replace a® = a\ = \/No ~ VN. In words, 
one assumes that the zero mode is macroscopically occupied and all commutators 
[ao, %] in the Hamiltonian are suppressed by relative powers of l/N; the quantum 
fluctuations of the zero mode may therefore be neglected. This, in combination with 
taking into account the constraint 

N = N + Y,a{a k (16) 

leads to the Hamiltonian 

U = Yl ( k2 ~ aN l 2 ) a ^ ak ~ \ aN Yl { a i a -k + a k a -k) + 0(1/N) . (17) 

All interaction terms are suppressed by l/N and go to zero in the double scaling 
limit iV — > oo, a — > with aN finite. The Hamiltonian can be diagonalized 

n = J2 ^b{b k , e k = y/k*(k 2 -aN) (18) 
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with a Bogoliubov transformation 

a k = u k b k + vlfr_ k , 
where the Bogoliubov coefficients are 



1 + 



L2 aV 
™ 2_ 

efc 



«2 



i 



qjV 
2 



(19) 



(20) 



The Boguliubov approximation breaks down whenever an e k becomes too small. In 
that case the initial assumption that only the zero mode is macroscopically occupied 
is no longer justified. Obviously, it is e\ that first goes to zero, namely when aN — > 1. 
Right at the phase transition, the Boguliubov approximation is never valid. It is 
worth noting however, that for any finite distance S from the critical point, there 
exists a minimal iV for which the approximation is valid. In other words, for any 
finite 5, the Boguliubov approximation becomes exact in the limit N — > oo. This is 
due to the fact that both the interaction terms as well as v\/N vanish in this limit 
for any finite 5. For 5 = 0, however, this is never true. 

In the aN > 1 case, the classical background is not homogenous any more, but 



is given by the bright soliton solution (12). In this case, the background induces 
an additional nontrivial mixing between momentum eigenmodes of different \k\. A 
decomposition into momentum eigenmodes requires an (unknown) analytic expres- 
sion for the Fourier components of the soliton and is thus no longer convenient. 
On the other hand, an analytic Bogoliubov treatment is still possible by directly 
decomposing Sip into normal modes: 



E 



u. 



If the mode functions obey the Boguliubov-de Gennes equations 



d e Uj + a^ (2uj + vj) + fj,Uj 
dgVj + a$>^{2vj +Uj) + fiUj 



E jUj 
-EjUj 



(21) 



(22) 
(23) 



and are normalized such that they form a complete set and the transformation (21) 



is canonical, the Hamiltonian is diagonalized. The first excited Bogoliubov modes 
have the form 



Ul (6) = iVisn 2 
Vl (9) = -iVicn 2 



The coefficient N% is defined by 



K{m) 

71 

K(m) 

7T 



mK(m) 



0n 



m 



m 



2tt [(2 - m)K(m) - 2E(m) 



(24) 
(25) 

(26) 
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2.3 Numerical Diagonalization 

While the Bogoliubov treatment provides an approximative description of the Bose 
gas deep in the respective phases, it fails, as we have reasoned above, around the 
critical point. 

A complementary method to explore the quantum properties of the system is 
numerical diagonalization of the Hamiltonian. Of course, numerical techniques are 
only applicable for sufficiently small Hilbert spaces. The Hamiltonian ^ is number 
conserving. This allows for exact diagonalization of ^ by considering a subspace of 
fixed N. However, to make any numerical procedure feasible, we need to limit the 
allowed momenta. In the spirit of [3], we truncate the basis of free states in which we 
perform the diagonalization to \l\ = 0, 1. This gives a very good approximation to 
the low energy spectrum of the theory well beyond the phase transition. Analytically, 



this can be seen by analyzing the spectrum of the soliton solution (12). Only for 
aN > 1.5, higher I modes start giving relevant contributions. We have further 
verified this numerically by allowing for \l\ = 2,3; the low energy modes are only 
marginally affected up until aN ~ 2. Our code allowed us to consider particle 
numbers N ^ 10000. In order to illustrate scaling properties, all analyses are 
performed for various particle numbers. 

Since the normalized coupling aN is the relevant quantity for a phase transi- 
tion, one can analyze all interesting properties for a fixed N by varying a. The 
corresponding spectrum of excitations above the ground state as a function of aN 
is shown in Fig |T] for iV = 5000 and — 1 < k < 1. One observes a decrease in the 
energy gap between the low lying excitations due to the attractive interactions as 
aN is increased. At the quantum critical point, the spacing between levels reaches 
its minimum. Its magnitude depends on the particle number N; the energy of the 
lowest lying excitation decreases with N. By further increasing the coupling a one 
reaches the solitonic phase. The spectrum corresponds to that of translations and 
deformations of a soliton. 

Obviously, ^ is invariant under translations; since we are considering a finite 
length ring, the ground state obtained by exact diagonalization can never correspond 
to a localized soliton. It will instead contain a superposition of solitons centered 
around arbitrary 9. This problem can be overcome by superposing a weak symmetry 
breaking potential to break the degeneracy between states with a different soliton 
position: 

H sb = H + V e (27) 

Ve = ^ J deft(e)co8dj>{e). (28) 

The higher e, the deeper the symmetry breaking potential, and the more localized 
the soliton will be. 
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0.0 0.5 1.0 1.5 2.0 

Figure 1: Energy spectrum for N = 5000 as a function of the effective coupling aN 



3 Quantum Phase Transition in the lD-Bose gas 

The mean field treatment of the attractive ID Bose gas above has signalled a quan- 
tum phase transition. The degeneration of the Bogoliubov modes at aN = 1 sup- 
ports the existence of a critical point. Although, by definition, a phase transition 
can only occur for infinite N, indications for it should already be visible for large 
but finite N. Here we will focus on two observations: 

(i) The one-particle entanglement entropy displays a sharp increase close to the 
critical point. 

(ii) The ground state fidelity peaks at the critical point; the height of the peak 
grows with N. 

3.1 One-Particle Entanglement 

The one particle entanglement entropy is defined as the von Neumann entropy Si = 
Tr[p lp log pi p ] of the one particle density matrix p lp of the ground state, obtained by 
singling out one particle and tracing over all iV — 1 other. As long as the ground state 
of the system is well described by a Hartree, i.e. product state, p± p describes a pure 
state; the entanglement entropy vanishes. When the critical point is approached, 
collective effects become important. No longer is the ground state described by a 
product state; consequently the entanglement entropy increases - a single particle 
becomes strongly entangled with the rest of the system. 
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Figure 2: One-particle entanglement entropy for N = 50 (blue, dotted), 500 (black, 
dashed), 5000 (red, solid). 



The one particle density matrix is defined via 



p lp = Tr p= Tr |0 GP )(0 GP |, (29) 

(JV-l)P (iV-l)p 



or, explicitly, in the one particle momentum eigenbasis 

N 

Here, is the occupation number of the k-th momentum mode and we have used 



{Pip)ij = ^j^2\a{n k }\ 2 ^- (30) 

{n k } 



\°gp) = ^2 a {n k }\{n k }) • (31) 
{«fc} 

We have plotted the numerically evaluated one particle entanglement as a function of 
aN for different N in Figj2j The increase close to the critical point gets profoundly 
sharper for larger N. Independent of N, the entropy is bounded by S ma ^ = log 3, 
due to the truncation of the one-particle Hilbert space to a three level system. 

The entanglement entropy becomes maximal for large aN. This, as argued 
before, is due to the fact that the numerical groundstate is given by a superposition 
of solitons localized at arbitrary positions]^] 



2 This quantum behavior is not expected to survive in the large N limit if the symmetry breaking 
potential is turned on. In this case, a vacuum is selected which does not mix with translated states 
(see Appendix). The entanglement entropy is therefore much smaller. 
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Figure 3: Numerical ground state fidelity susceptibility for iV = 3000 (blue, dashed), 
N = 5000 (black, dotted) and N = 10000 (red, solid). 



3.2 Ground State Fidelity 

Ground State Fidelity (GSF) was introduced in [7] as a characteristic of a QPT. It 
is defined as the modulus of the overlap of the exact ground states for infinitesimally 
different effective couplings. 

F(aN,aN + 5) = |(0 aJV |0 aJV+5 )| (32) 

Far away from the critical point, this overlap will be very close to unity. For small 
aN, the ground state is dominated by the homogeneous state, and while coefficients 
may change slightly, no important effect will be seen. The analogous statement 
holds deep in the solitonic regime. While the shape of the soliton changes, it will 
so smoothly; in the infinitesimal limit, the overlap is one. Right at the critical 
point, however, the ground state changes in a non-analytic way. The homogeneous 
state ceases to be the ground state and becomes an excited state, while the soliton 
becomes the new ground state. As energy eigenstates with different eigenvalue are 
orthogonal, the ground state fidelity across the phase transition is exactly zero. 

The GSF has the disadvantage of depending on the arbitrary choice of the small 
parameter 5. This can be cured by introducing the fidelity susceptibility Xg S (aiV) 
as the second derivative of the GSF. 

F(aN,aN + 5)-F(aN,aN-5) 
Xgs (« AO = hm . (33) 

It has been shown [8] that singular behavior of the fidelity susceptibility directly 
signals a discontinuity of the first or second derivative of the ground state energy - 
a quantum phase transition. 
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The aforementioned behavior is of course idealized for an infinite system, where 
ground state degeneracy and thus level crossing become an exact property In the 
finite iV systems we examined numerically, the overlap cannot go to zero, because 
there is anticrossing which allows the energy levels to degenerate only for iV — > oo. 

Still we can observe a drop in the fidelity which deepens with N but is of magni- 
tude much smaller than 1 for all N we were able to simulate. The fidelidy suscepti- 
bility as obtained from the exact diagonalization is plotted in Fig|3]for different N. 
In the limit iV — » oo, we expect a behavior x gs — > —5"(aN — 1). This tendency can 
be clearly observed. Both the negative and the positive peak move towards aN = 1, 
they become narrower, and their modulus diverges with growing N. 



4 Fluctuation Entanglement 

We will now consider the entanglement between the fluctuation Sdk = dk — a c k of a 
given original momentum mode and the fluctuations of the rest of the system. The 
motivation for studying this quantity is twofold. We imagine, that an external ob- 
server would couple linearly to the bosonic field (so that the situation has some min- 
imal resamblance with the gravity case). It has been pointed out [9] that for such a 
coupling, field values (or their Fourier components) will be the environment-selected 
pointer states^ and not localized single particle states. This leads us to consider the 
entanglement of a momentum mode, rather than single-particle entanglement, as 
a measure of relevant quantum correlations of the given state. Furthermore, the 
observer couples to the original field dk and hence its fluctuations as opposed to 
coupling to the Bogoliubov modes b\. 

More technically speaking, the quantity we calculate is the von Neumann entropy 
of the reduced density matrix for a given Sdk 

(Sp k )nm= Tr \p(6alr\0 c )(0 c \(6a k r} (34) 

modes k'^pk L 

where |0 C ) denotes the state that would be observed classically. 

Fluctuation entanglement provides a measure for the quantum correlations be- 
tween a single momentum mode with the rest of the system. It hence gives a direct 
handle of the quantumness of our ground state as measured by an outside observer 
if coupled linearly to the field. Note also that due to the fact that we are considering 
a closed system, the fluctuation entanglement is exactly equivalent to the Quantum 
Discord introduced in works \±0\ as a measure of quantumness. 



3 Pointer states denote those states that are stable with respect to interactions with the envi- 
ronment and therefore correspond to classically observable states. 
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4.1 Calculation in the Bogoliubov Approximation 



In order to calculate the fluctuation entanglement in the Bogoliubov case, note that 
the sought-for density matrix is GaussiariQ The ground state in terms of b k is 
Gaussian and the Bogoliubov transformation amounts to squeezing - which leaves a 
Gaussian state Gaussian. Also integrating out modes in a Gaussian state does not 
change this property. Hence the reduced density matrix in terms of 8a k must have 
the form 



Pk = C k exp 



-A fe ( 8a\8a k + -^r k 



1 

2 



8a) k 8a} k + 8a k 8a k 



with real coefficients A& and and normalization C k such that Trp^ 
density matrix has a von Neumann entropy 



(35) 



1. This 



Sk 



Tt 



coth 



A fcv /r 



In 1 



(36) 



We can fix the unknown coefficients by imposing 



(i/j\Sa[Sd k \ip) = Tr[p k 5al5a k ] 



and 



(ip\8a k 8a k \ip) = Tr\p k 8a k 8a k ], (37) 



where is the groundstate of the Bogoliubov modes. 



4.2 Homogenous Phase 



In the homogenous case, imposing (37) and evaluating the left hand side by inserting 
the Bogoliubov transformation ( 19 ) leads to 

A fc = ln^^ , r k = and C k = l/u 2 k . (38) 

Thus, the fluctuation entanglement entropy is 

S k = u 2 k \nul-v 2 k \nv 2 k . (39) 

The entanglement of the first momentum mode Si diverges near the critical point 

aN = 1 — 8 as 

Si » l-ln(4) - i]n& (40) 

A similar divergence of an entanglement entropy has been pointed out in spin chain 
(and analogous) systems undergoing a phase transition [TT]. In contrast to these 
cases however, where the entanglement is between nearest neighbour sites, the di- 
verging entanglement in our case is between different low-momentum modes and 



4 A density matrix is called Gaussian, when its Wigner function W(a,a*) 
^2 / d 2 /3exp(—if3a* — i(i*a)Tv[pexp(i(3a 1 ' +i(3*a)] is Gaussian 
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Figure 4: Analytical fluctuation entanglement 



not between localized sites. So one may say, that the entanglement in our case is 
long-range. Furthermore it should be noted, that the entanglement of the higher 
modes \k\ > 1 stays finite near the critical point, showing that the diverging en- 
tanglement is an infrared effect, which can be expected to be independent of short 
distance physics. 

4.3 Solitonic Phase 

The relevant expectation values in the Bogoliubov ground state are given by 

(*lj\5ai5aM = E (/ ^V^) (/ e- irU, v k (6)*(W) , (41) 

^\5a m 5aM = E (/ ^ m9 u k {e)d6^ (^j e~ int 'v k {6y 'd$j . (42) 

It can be checked that close to the phase transition the first excited mode gives the 
leading contribution to the aforementioned entanglement entropy. The quantities 
{il)\5a\5a\\ip) and {^\8a\8a\\il)) can be obtained by numerical integration. The pa- 
rameters A, t of the reduced gaussian density matrix can then be determined. The 
final von Neumann entropy again shows a divergence^] close to the phase transition. 
Fig|4] shows the fluctuation entanglement obtained in the Bogoliubov approximation 
on both sides of the phase transition. 



5 Thc shape of the divergence obtained by numeric integration seems to be consistent with a 
logarithm with a coefficient close to 0.33 near the phase transition. 
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4.4 Numerical Treatment 

As discussed before, for any given finite size system, the Bogoliubov approximation 
should not be trusted close to the critical point. Therefore it is important to study 
the exact behavior of finite size systems numerically, in order to substantiate the 
claim that the fluctuation entanglement entropy becomes large. 

Within an exact treatment this quantity is considerably more difficult to extract, 
because in contrast to the Bogoliubiv analysis, one does not have direct access to 
a "classical background" which one could use to disentangle classical correlations. 
Instead, the seperation can be obtained through the following procedure. 

Since all numerical solutions are obtained for a fixed particle number N, the field 
expectation value in the exact ground state |0jy) will necessarily vanish, (0 N \ip(9) |0jv) : 
0. Obviously, the ground state |0jy) can hence never correspond to the classical 
(coherent) state with a wave function corresponding to the soliton solution of the 
Gross-Pitaevskii equation, (V'ci|'0(^)|'0ci) = ^gp($)- 111 order to define a mapping 
from |0jv) to iV'ci)) we numerically search for the coherent state |a) with maximal 
overlap with \0n). This state is expected to be annihilated by the perturbations of 
the Gross-Pitaevskii ground state, 

5a k \a k ) = 0, (43) 

where 

Sa k = a k - c k (44) 



and the c k are the Fourier coefficient of ^/gp(9). From Eq.(43) it directly follows 



that a k = c k . There is now an obvious measure of correlations which excludes those 
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of the Gross-Pitaevskii background: The entanglement entropy of the 8a\ modes, 
described by the density matrix 



l)kl 



Tr [p|5Zi)(5fci|] 



(45) 



Now, \5l\) denotes the eigenstate of 8N\ = 5a[5ai with eigenvalue 5l\. Eq. (45) 



directly corresponds to the density matrix (35) in the Bogoliubov approximation 



Using the relations (43) and (44), it can be directly recast to take on the form 



(Pi)jh = Tr 



|ai)(ai| 



\ai - ax) 



(46) 



which, by the definition of a coherent state, can be straightforwardly evaluated. 

The resulting fluctuation entanglement is shown in Fig|5] for different particle 
number. It has a clear maximum at the would-be-phase-transition. The maximum 
value becomes larger and the peak narrower with increasing particle number, so 
the divergence in the Bogoliubov case seems a plausible limit. The fact that at 
aN = 2, the fluctuation entanglement is still quite high is not surprising. Only in 
the limit iV — > oo do we expect to see the behavior of the Boguliubov analysis. This 
is supported by the fact that a decrease is observed for increasing N, as well as for 
stronger localization potentials. 



5 Conclusions and Outlook 

In this paper, we have considered properties of the 1+1-dimensional attractive Bose 
gas around its critical point. By analyzing important indicators for QPTs, we pro- 
vided further evidence that a tuning of the effective coupling gN leads to a phase 
transition in the system. More importantly, we have shown that quantum corre- 
lations become very important close to the critical point - contrary to the naive 
intuition that at sufficiently large particle number, systems should behave approxi- 
mately classical. We have also pointed out that the quantum entanglement of the 
bosons close to the critical point is "long range" - in contrast to the observations in 
spin-chain systems that display nearest neighbour entanglement at criticality. 

The motivation for our study of this model system, however, was the conjecture 
that black holes are bound states of a large number of weakly interacting gravitons. 
It has been claimed that the graviton condensates behave significantly different with 
respect to the semiclassical black hole analysis due to their being at a quantum criti- 
cal point. It was argued that criticality allows quantum effects to only be suppressed 
by the perturbative coupling a g ~ 1/N. If the qualitative insights from our simple 
toy model are valid for graviton condensates our results can back up several of the 
claims. We can argue that quantum effects become important for attractive Bose 
condensates at their critical point - even though the perturbative coupling is very 
small. And that the entanglement of the true state is long range - consistent with the 
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notion of a condensate of gravitons of wavelength comparable to the Schwarzschild 
radius. This would imply that for a black hole, the semiclassical treatment with 
a background geometry that obeys classical general relativity and quantization of 
fields on top of this rigid background becomes invalid much earlier than what the 
standard lore tells. Although curvature invariants in the horizon region of large 
Schwarzschild black hole are small, the semiclassical treatment is not applicable. 
Instead, quantum correlations in the graviton bound state become relevant. Im- 
portantly, our results point in the direction that the physics is dominated by large 
wavelengths. Therefore the description of black holes as graviton condensates has 
the attractive feature of being independent of the ultraviolet completion of gravity. 
The only requirement being that the low energy theory resembles perturbatively 
quantized Einstein theory with a massless spin two graviton. 

The 1+1-dimensional Bose gas can indeed capture quite a few of the intriguing 
features of black holes and their possibly quantum nature. To understand in more 
detail time dependent features, such as Hawking evaporation, resolutions of the in- 
formation paradox or scrambling, the implementation of dynamical methods will be 
amongst the aims of immediate future work. This will then necessarily also address 
possible couplings to external systems in order to be able to model the evapora- 
tion process. Working with more spatial dimensions may prove feasible to model 
the collapse induced by Hawking evaporation. This could alternatively be achieved 
by considering couplings that show further resemblance with gravitational self- in- 
teractions. Steps in these directions also include generalizations to non-number 
conserving, and ultimately relativistic theories. 

Instabilities can in turn be countered by adding repulsive interactions that dom- 
inate at very short scales. Stable configurations of that sort would correspond to 
extremal black holes. Their properties also provide a vast playground for future 
investigation. 
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A Appendix 

A.l Spontaneous symmetry breaking in finite volume 

Standard lore states that there can be no spontaneous symmetry breaking in finite 
volume with a finite number of fields. This may seem puzzling since we claim here 
explicitly that our ground state is a localized soliton for gN > 1, clearly signaling 
spontaneous breaking of translation invariance. This puzzle is resolved by noticing 
that the mean field approximation can only be exact in the limit N — > oo. 

To illustrate this, consider a number conserving Hamiltonian containing a finite 
number 2m of fields in the potential. Mean field approximation becomes exact if the 
corresponding state can be approximated arbitrarily well by a Hartree-type state, 
i.e. 

^ ) = ®J2c k \k). (47) 

k 

A translated state is given by 

k 
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Now it can be seen that spontaneous symmetry breaking is indeed possible as we 
take the limit N — > oo. It can be shown that 



2N 




(49) 



k 



2(N-m) 




(50) 



k 



In both expressions the right hand side vanishes as long as ^ for at least 
two different k. In other words, in the large N limit localized objects, localized at 
different points are orthogonal and, furthermore, do not mix under time evolution. 

In the proof for absence of finite volume symmetry breaking enters the assump- 
tion that expectation values of composite operators made out of the fields are finite. 
This assumption breaks down in the large N limit. While this implies unbounded- 
ness of the energy in this limit, it is nothing to worry about. The energy diverges 
linearly in N and hence the energy per particle remains bounded, reminiscent of a 
thermodynamic limit. 
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